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o : 

' Abstract- We find a new class of exact solutions to the Einstein-Maxwell 

i— 1 1 equations which can be used to model the interior of charged relativistic ob- 

Qi jects. These solutions can be written in terms of special functions in general; 

for particular parameter values it is possible to find solutions in terms of 
'— elementary functions. Our results contain models found previously for un- 

i ■ charged neutron stars and charged isotropic spheres. 

O ' Keywords- charged spheres, Einstein-Maxwell equations, relativistic astro- 

physics. 

00 : 1- INTRODUCTION 

o : 

The Einstein-Maxwell system of field equations are applicable in mod- 
elling relativistic astrophysical systems. We need to generate exact solutions 
to these field equations to model the interior of a charged relativistic star 
that should be matched to the Reissner-Nordstrom exterior spacetime at the 
boundary. A general treatment of nonstatic spherically symmetric solutions 
with vanishing shear was performed by Wafo Soh and Mahomed pQ using 
symmetry methods. The matching of nonstatic charged perfect spheres to 
the Reissner-Nordstrom exterior was considered by Mahomed et al. [2] who 
showed that the Bianchi identities restrict the number of solutions. Partic- 
ular models generated can be used to model the interior of neutron stars as 
demonstrated by Tikekar [3j, Maharaj and Leach jl] and Komathiraj and 
Maharaj [5]. Charged spheroidal stars have been widely studied by Sharma 
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et al. [B] and Gupta and Kumar |7j. There exist comprehensive studies 
of cold compact objects by Sharma et al. [8], analysis of strange matter 
and binary pulsar by Sharma and Mukherjee [9] and quark-diquark mix- 
tures in equilibrium by Sharma and Mukherjee [ID] , in the presence of the 
electromagnetic field. Charged relativistic matter is important in the mod- 
elling of core-envelope stellar systems as demonstrated by Thomas et al. [TTj . 
Tikekar and Thomas [12] and Paul and Tikekar [13] . The recent treatment of 
Thirukkanesh and Maharaj [H] deals with charged anisotropic matter with 
a barotropic equation of state which is consistent with dark energy stars and 
charged quark matter distributions. 

The exact solution of Tikekar [3] is spheroidal in that the geometry of the 
spacelike hypersurfaces generated by t=constant are that of a 3-spheroid. 
This condition of a spheroid helps to mathematically interpret the solu- 
tion since it provides a transparent geometrical interpretation. On physi- 
cal grounds we find that this solution can be applied to model superdense 
stars with densities of the order 10 14 g cm 3 . The physical features of the 
Tikekar model are therefore consistent with observation, and consequently 
it attracts the attention of several researches as a realistic description of the 
stellar interior of dense objects. This solution was extended by Komathiraj 
and Maharaj [5] to include the electromagnetic field, with desirable phys- 
ical features. In this paper we show that a wider class of solutions to the 
Einstein-Maxwell system is possible by adapting the form of the gravitational 
potentials. Our intention is to obtain simple forms for the solutions that are 
physically reasonable and may be used to model a charged relativistic sphere. 

2. SPHERICALLY SYMMETRIC SPACETIME 

The metric of static spherically symmetric spacetimes in curvature coor- 
dinates can be written as 

ds * = _ e 2u(r) dt 2 + e 2X(r) dr 2 + ^2 + ^2 ^ ^ 

where z/(r) and A(r) are two arbitrary functions. For charged perfect fluids 
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the Einstein-Maxwell system of field equations is given by 

1 /. _o\x 2 A' _ \ I 
r 2 v r 2 



(1 - e- 2A ) + — e- 2X = p + -E 2 (2a) 
r 2 

4(l-e" 2A ) + — e- 2A = P-^ 2 (2b) 

e~ 2X (v" + v' 2 + --v'\'--) = p + -E 2 (2c) 
\ r r J 2 

a = ^e-\r 2 E)' (2d) 

for the line element (pQ). The quantity p is the energy density, p is the 
pressure, E is the electric field intensity and a is the proper charge density. 
To integrate the system ([2} it is necessary to choose two of the variables 
u, A, p, p or 2£. In our approach we specify A and E. 
In the integration procedure, we make the choice 

2A ^ / /o\ 

where and / are arbitrary constants. Note that the choice ([3]) ensures 
that the metric function e 2X is regular and finite at the centre of the sphere. 
When k = —7 and I — 1, in the absence of charge, we regain the Tikekar 
interior metric [3] which models a superdense neutron star. Also observe 
that when / = 1 we regain the metric function considered by Komathiraj 
and Maharaj [5] which generalises the Maharaj and Leach [1] and Tikekar 
[3] models. Therefore particular choices of the parameters k and I produce 
regular charged spheres which are physically reasonable. Also the choice ([3]) 
ensures that charged spheres generated, as exact solutions to the Einstein- 
Maxwell system, contain well behaved uncharged models when E = 0. On 
eliminating p from (12bl) and fl2c|) , for the choice ([3]) , the condition of pressure 
isotropy with a nonzero electric field becomes 

(l-A;r 2 AR 2 ) 2 £ 2 = (l - kr 2 /R 2 ) (l - lr 2 / R 2 ) fu" + v' 2 - ^ 

-(/-fc)^(,' + i) + ^(l-^ 2 ) (4) 
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which is nonlinear. 

To linearise the above equation it is now convenient to introduce the 
transformation 

ip(x) = e u (5a) 
x 2 = l-lr 2 /R 2 (5b) 

where 1^0. This transformation helps to simplify the integration proce- 
dure but changes the form of the potentials and matter variables. Then (Tj0) 
becomes 

2 , •• . , f(l-k + kx 2 ) 2 R 2 E 2 k(k-l) 



(l-k + kx l )i) - kxip + 2 J + — J -J i) = (6) 

in terms of the new dependent and independent variables ip and x respec- 
tively. Equation dSJ) must be integrated to find if), i.e. the metric function A. 
Note that the Einstein-Maxwell system (j2j) implies 

l(l-k)(3l-k + kx 2 ) 1 a 

p = __ \ — — ( -2lx— +k — l\ + -E 2 (7b) 

F R 2 (l-k + kx 2 ) I V / 2 

= / 2 [2^-(l-xW 

i? 2 (l-:r 2 )(Z-fc + /b; 2 ) 1 J 

in terms of the variable x. Note that we have essentially reduced the solution 
of the field equations to integrating (jHJ). It is necessary to specify the electric 
field intensity E to complete the integration. Only a few choices for E are 
physically reasonable and generate closed form solutions. We can reduce © 
to simpler form if we let 

2 = akl{x 2 - 1) 

R 2 (l-k + kx 2 ) 2 { ' 

where a is constant. When a = or A; = there is no charge. The form 
for E 2 in (IS]) vanishes at the centre of the star, and remain continuous and 
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bounded in the interior of the star for a wide range of values of the parameters 
a, k and I. Upon substituting the choice (|HJ) into (jSJ), we obtain 



which is the master equation for the system (J7)). We expect that our investi- 
gation of equation (jUj) will produce viable models of charged stars since the 
special cases a = and cr ^ 0, k ^ 0, 1 = 1 yields models consistent with 
neutron stars. 



As the point x = is a regular point of (JUJ), there exists two linearly 
independent series solutions with centre x = 0. Thus we must have 



where aj are the coefficients of the series. For an acceptable solution we need 
to find the coefficients a,i explicitly. On substituting (|10|) in (J2J) we obtain 
after simplification 

1(1 - k)(i + l)(z + 2)a i+2 + k[a + k-l + li(i - 2)]^ = 0, i > (11) 

The equation (ITT)) is the basic difference equation governing the structure of 
the solution. It is possible to express the general form for the even coefficients 
and odd coefficients in terms of the leading coefficient a$ and a\ respectively 
by using the principle of mathematical induction. We generate a pattern 



1(1 — k + kx 2 )4> - klxtp +k(k-l + a)ip = 



(9) 



3. NEW SOLUTIONS 




(10) 



i=0 




L_ Y[[a + k-l + l(2q - 2)(2q - 4)]a (12) 



for the even coefficients ao, a,2, 04 . . .. Also we find the pattern 



a 2i+l — 




) fib + A; - Z + Z(2g - l)(2g - 3)] ai (13) 



5 



Mathematical and Computational Applications 
©Association for Scientific Research 

for the odd coefficients a±, a 3 , . . .. Here the symbol J7 denotes multipli- 
cation. 

From ( TTUj) . (|T2|) and (|T3|) . we can write the general solution of <Q as 



Thus we have found the general solution to the differential equation ([9]) for 
the particular choice of the electric field ((HJ). Series (I15a.[) and ( 115bl) converge 
if there exists a radius of convergence which is not less than the distance 
from the centre of the series to the nearest root of the leading coefficient in 
(J9j). This is possible for a range of values of k and I. 

The general solution ( fl4l) is given in the form of a series which may be 
used to define new special functions. For particular values of the parameters 
a, k and I it is possible for the general solution to be written in terms of 
elementary functions which is a more desirable form for the physical descrip- 
tion of a charged relativistic star. Solutions that can be written in terms 
of polynomials and algebraic functions can be found. This is a lengthy and 
tedious process and we therefore do not provide the details; the procedure 
is similar to that presented in Komathiraj and Maharaj [5] which can be 
referred to. The solutions found can also be verified with the help of soft- 
ware packages such as Mathematica. Consequently we present only the final 
solutions avoiding unnecessary details. 




(14) 



where we have set 




(15b) 



(15a) 
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Two classes of solutions in terms of elementary functions can be found. 
These can be written in terms of polynomials and algebraic functions. The 
first category of solution for if)(x) is given by 



n-2 



{n-j -2)!(2j + l)! 

with the values 

7 = 4 

4ln(n - 1) + a 

k + a = /[2-(2n-l) 2 ] 
The second category of solution for ip(x) has the form 

«*) - A EH ^j!ti".' l ' v ^' 



(n-j)!(2i + l)!' 



n-1 



- fc + A;* 2 ) 3 / 2 E(-^) % (W+ iL/ ( 17 ) 

(n-j - l)!(2j)! 

with the values 

, 4/ 

^ 4/n 2 - / + a 

k + a = 2/(1 -2n 2 ) 

where A and -B are arbitrary constants and x 2 = 1 — lr 2 /R 2 . 

5. SPECIAL CASES 

From our general class of solutions ( 1T61) and (fTTI) . it is possible to generate 
particular cases found previously. These can be explicitly regained directly 
from the general series solution (fl4"l) or the elementary functions (Tl6|) and 
(jT7|) . We demonstrate that this is possible in the following classes. 
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We set k + a = —7l(n = 2). Then 7 = 4(7/ + a)/ (81 + a) and it is easy 
to verify that equation (!T6|) becomes 



71 -r (\ \ 2 8 (11 + a\ 2 4 



^(x) = A' 1 - 4 — x 2 + . , 



+E'x 1 



7/ + a\ 2 ^ 3/2 



81 + a 



x 



where A' = A/2 and 5' = 25(8/ + af' 2 are new constants. Further setting; 
a = and / = 1, we obtain 

7 „ 49 A „, / 7 ^ 3/2 



= A' 11 - -x 2 + -x 4 J + S'x ^1 - -x 2 J (18) 

and x 2 = 1 — r 2 /R 2 . Thus we have regained the Tikekar model [3 J which is 
a viable model in the modelling of superdense stars. 

We set k + a = -2l(n = 1). Then y, = 4(2/ + a)/ (31 + a) and (E 
becomes 



+ B" l 1 /2( + QN ^ 



3/ + a 



x 2 



where are A" = A and B" = B(3l + a) 3 / 2 are new constants. Further setting 
a = and l/R 2 = C/2 (k/R 2 = -C) and letting X = Cr 2 we obtain 

i, = ^=(2 - X)^( 5 + 2X) + J^(l + X) 3 / 2 

where we have set ip = ip(X). Thus we have regained the Durgapal and 
Bannerji [15] model which is widely used in the modelling of neutron stars. 

If we set / = 1 and a = then f|T6|) and (IT7|) reduce to the corresponding 
expressions in the solution of Maharaj and Leach [4] which implies a wide 
family of models for uncharged relativistic spheres which have the advantage 
of being expressed in elementary functions. 
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If we set I = 1 then (fTSj) and f[T7|) contain the solution of Komathiraj and 
Maharaj [5] for charged spheres which are generalizations of earlier models 
with spheroidal geometry. 

5. DISCUSSION 

We have studied the Einstein-Maxwell system of equations for a particular 
choice of the electric field intensity. The gravitational potential was gener- 
alised to include the spheroidal geometry of the hypersurfaces t=constant of 
previous investigations. When I = 1 then we regain the Tikekar [5] model 
and other exact solutions found previously. We demonstrated that it was 
then possible to reduce the condition of pressure isotropy to a second order 
linear ordinary differential equation. This equation can be solved in general 
using the method of Frobenius and the solution are in terms of new special 
functions. Solutions in terms of elementary functions can be extracted from 
the general solution for specific parameter values. Particular models studied 
previously are contained in our general class of solution. These solutions 
may be useful in studying the physical behaviour of dense charged objects in 
relativity which will be the objective in future work. 

We briefly discuss the behaviour of the matter variables close to the cen- 
tre. We can graphically represent the the matter variables in the stellar 
interior for particular choices of the parameter values. To this end we have 
produced Figure [1] with the help of the software package Mathematica. We 
have set A — B — C — 1, k — — |, I — — 1 and a = | over the interval 
< r < 1, to generate the relevant plots in Figure [TJ Plots A and B denote 
the profiles of energy density p and the pressure p; plot C denotes the elec- 
tric field intensity E 2 . We observe that these matter variables remain regular 
in the interior. We note that the energy density p and the pressure p are 
positive and finite; they are monotonically decreasing functions in the inte- 
rior. The electric field intensity E 2 is positive and monotonically increasing 
in this interval. Thus the quantities p, p and E 2 are finite, continuous in the 
interval. 
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Figure 1: Plots of the matter variables p, p and E 2 . 
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